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Objects that float at the interface between a liquid and a gas interact because of interfacial 
deformation and the effect of gravity. We highlight the crucial role of buoyancy in this interaction, 
which, for small particles, prevails over the capillary suction that is often assumed to be the dominant 
effect. We emphasize this point using a simple classroom demonstration, and then derive the physical 
conditions leading to mutual attraction or repulsion. We also quantify the force of interaction in 
some particular instances and present a simple dynamical model of this interaction. The results 
obtained from this model are then validated by comparison to experimental results for the mutual 
attraction of two identical spherical particles. We conclude by looking at some of the applications 
of the effect that can be found in the natural and manmade worlds. 



I. INTRODUCTION 

Bubbles trapped at the interface between a liquid and 
a gas rarely rest. Over a timescale of several seconds to 
minutes, long-lived bubbles move towards one another 
and, when contained, tend to drift towards the exterior 
walls - as shown in fig.[l] The sceptical reader may read- 
ily verify these claims by pouring themselves a glass of 
sparkling water (or, if they prefer, wine) and following 
the motion of those bubbles at the surface - particularly 
those at the periphery of the glass. This phenomenon has 
even been affectionately dubbed the 'Cheerios effect' af- 
ter the observation that breakfast cereals floating in milk 
often clump together or stick to the walls of the breakfast 
bowl^ . 

In this article, we bring together a number of perspec- 
tives on the Cheerios effect gathered from the literature 
and our own experience at home, in the kitchen, and 
in the laboratory. We show how simple physical ideas 
lead to an understanding not only of the attraction it- 
self, but also of its dynamical consequences. Despite be- 
ing a subject with enormous potential for simple, reliable 
party tricks, the technological implications of the Chee- 
rios effect are far from frivolous. Indeed, much research is 
currently being undertaken to investigate the possibility 




FIG. 1: Bubbles floating on water in a petri dish. The bubbles 
are observed to aggregate before moving to the wall of the 
container. After sufficiently long times, the island of bubbles 
in the centre also migrates to the wall. 



of using surface tension to induce the 'self-assembly' of 
small-scale structures'^. Understanding the way in which 
particles aggregate at an interface and thence being able 
to control the form of the aggregate as well as the dynam- 
ics of its formation may one day enable much simplified 
manufacture of components of micro-electromechanical 
systems (MEMS). 

For floating objects in equilibrium or motion, we must 
consider the balance of linear momentum both in the 
flotation plane as well as out of the plane, and in addition 
the balance of angular momentum in all three directions. 
Many of the misconceptions in the field arise from consid- 
ering only some but not all of these balance equations. 
In particular, neglecting the condition of vertical force 
balance leads to an underestimation of the importance of 
the particle's buoyancy in determining the nature of the 
interaction. 

We begin with a discussion of the physical mechanism 
that leads to the observed attraction in most instances 
and illustrate the role of particle buoyancy by means of 
a simple experiment. We then focus on a series of sim- 
ple examples that allow us to quantify the magnitude of 
the attractive force. The first of these is inspired by an 
oversimplified physical picture that is often perceived as 
the dominant effect. By considering carefully the vertical 
force balance that must be satisfied for particles to float, 
we then show that it is this that is dominant for small 
particles and subsequently propose a simple dynamical 
model for the attraction of two spherical particles. Fi- 
nally, we show that consideration of the remaining equi- 
librium condition, that of torque balance, can lead to 
amphiphilic strips and conclude with some possible bio- 
logical implications of our observations. 



II. THE PHYSICAL ORIGIN OF ATTRACTION 

The mechanism behind the apparent attraction be- 
tween bubbles or between a bubble and the wall of a glass 
is easy to understand by considering the geometry of the 
interface at which the bubbles are trapped. For simplic- 
ity, we consider the latter case (schematically illustrated 
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FIG. 2: Schematic of a single bubble close to a wall, along 
with the definition of the contact angle. 

in fig. [2]), although the explanation of the clustering of 
many bubbles is similar. Here, the air-water interface is 
significantly distorted by the presence of the wall (the 
well-known meniscus effect) and, since the bubble is pos- 
itively buoyant, there is a net upward force, Fg, on the 
bubble. Because it is constrained to lie at the interface, 
however, the bubble cannot simply rise vertically and in- 
stead does the next best thing by moving upwards along 
the meniscus. As glass is in general wetting (the contact 
angle 9, defined in fig. [H satisfies 9 < tt/2), in moving 
upwards along the meniscus the bubble also moves closer 
to the wall. Viewed from above, it appears as if there 
is an exotic attractive force acting between the wall and 
the bubble when in fact the bubble is merely obeying 
gravity and moving in a constrained path imposed by 
the presence of the wall. 

Now, a single bubble will deform the interface just as 
the presence of a wall does, though this time for a differ- 
ent reason and to a lesser extent. In the bubble's case, 
it can only remain at the interface because the buoy- 
ancy force (that is tending to push the bubble out of the 
liquid) is counterbalanced by the surface tension force 
(that opposes the deformation of the interface and hence 
is tending to keep the bubble in the liquid). These two 
competing effects reach a compromise where the bubble is 
partially out of the liquid but the interface is slightly de- 
formed. This deformation is significant enough to influ- 
ence other bubbles nearby, which move 'upwards' along 
the meniscus, and so spontaneously aggregate. 

This mechanism was first proposed by Nicolson^ as a 
means by which the bubbles that constitute a bubble 
raft interact and give the raft its solid-like properties. 
As we shall see in section |Vl this mechanism provides 
the dominant contribution for the interaction between 
sufficiently small particles. 

The same argument works in the case of particles that 
are significantly more dense than water floating at an 
interface. That such heavy particles can float at all is 
due to the fact that this time surface tension stops the 
interface from deforming too much downwards as would 
happen if the particle were to sink. This reversal of the 
interfacial curvature can be seen clearly in fig. [3] for a 
metal drawing pin floating on water: surface tension must 
act upwards to counterbalance the weight of the pin. By 
analogy with what is observed with bubbles, we should 




FIG. 3: A photograph of a drawing pin fioating upturned on 
water. Notice that the deformation of the interface in this 
case is the opposite to that around a bubble or near a wetting 
wall. 



expect that another drawing pin floating near the first 
will 'fair down the interface and hence the two appear to 
be attracted to one another - as is indeed observed. 



III. REPULSION: OFTEN MISUNDERSTOOD 

So far we have seen only that the deformation of an 
interface caused by the presence of particles at that in- 
terface can lead to mutual attraction between those par- 
ticles and eventually to the formation of large clusters. 
This, however, is only half of the story. Imagine that we 
were to fioat a buoyant bubble in the vicinity of a draw- 
ing pin - should they also attract? On the basis of the 
previous physical argument, we expect that the bubble 
will move upwards along the interface distorted by the 
presence of the drawing pin. However, since the pin is 
not buoyant (i.e. the interface has the curvature shown in 
fig. [3]), moving along the interface will, in this case, cause 
the bubble to move away from the drawing pin and so the 
two objects in fact repel one another. 

A more striking demonstration of this repulsion can be 
achieved using two drawing pins, provided they are the 
kind that has a thin plastic cap around the blunt end. 
As we should expect from the discussion given in section 
im these two drawing pins will attract when floated at 
the interface. However, if we now carefully remove the 
plastic cap from the top of one and float it (the cap) near 
the intact drawing pin, then the two will be seen to move 
apart. 

This simple experiment apparently challenges the com- 
mon assumption that the attraction or repulsion of par- 
ticles at interfaces should depend solely on the wetting 
properties, namely the contact angles, of the particles 
(see, for example, p. 70, ex. 3 of Batchelor's book" or 
Campbell et al^). Here, the wetting properties of the 
plastic cap (which is the only part that is in contact with 
the liquid) are not altered by removing it from the pin 
but the weight that it must support is much reduced, 
making the cap buoyant. In turn, this alters the balance 
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FIG. 4: The geometry of two infinite plates in a semi-infinite 
fluid. The planes have contact angles 61 and O2 respectively 
and are a horizontal distance d apart. 



between surface tension and gravity so that the interface 
must now pull down on the cap to keep it at the inter- 
face and so the deformation near the cap resembles that 
around a bubble. This change in the sign of the curvature 
of the interface has been brought about without chang- 
ing the surface properties of the cap, simply because of 
a change in the effective density of the particle: a possi- 
bility with many potential applications that appear not 
to have been explored fully, although it is well known to 
those who know well^. 



IV. 



A MODEL CALCULATION 



Quantifying the physical picture outlined in section HIl 
allows us not only to predict the conditions under which 
the interfacial curvature changes sign, but also to un- 
derstand at a simple level the dynamical interaction be- 
tween two particles. We start with an idealised problem 
in which we account only for the condition of horizontal 
force balance (and neglect the vertical force and torque 
balance conditions) by focusing on two infinite vertical 
plates at a liquid-gas interface, as shown in fig. [4l The 
presence of the plates distorts the interface leading to a 
force of attraction between the two plates whose mag- 
nitude we shall now calculate. This set-up has, in fact, 
been used as a model for explaining the Cheerios effect^ 
and while we argue later that this picture is incorrect 
for floating objects, it does lend itself to a simple cal- 
culation (for a further simplification leading to similar 
conclusions, the reader is referred to the Encyclopedia 
Brittanica article by Clerk-Maxwell^). 

The equation of the interface z — h{x) is determined 
from the condition that the pressure change across the 
interface due to surface tension (which is proportional to 
the curvature of the interface) is equal to the hydrostatic 
pressure difference caused by the deformation of the in- 



terface (see p. 65 of Batchelor-4 for a thorough discussion 
of this). For small interfacial deflections, this balance 
may be written as: 



(1) 



where 7 is the surface tension coefficient of the liquid-gas 
interface, p is the density of the liquid and g is the accel- 
eration due to gravity. Equation ([T|) is to be solved with 
the boundary conditions that the contact angles 9i and 
6*2 are given at each of the plates and that the deflection 
of the interface should decay far away from the plates. 
In each of the regions labelled i = 1,2,3 in fig. IH the 
solution of ll]) gives the equation of the interface to be 
z = hi{x) = AiC^^/^'' + BiC^/^" where = 99 is 
the capillary lengthy which defines the length-scale over 
which interactions occur. The conditions /ii(— 00) = = 
h^{oo) give Ai = = -83, which, combined with the con- 
tact angle conditions h[ (0) — cot 9i and h'^{d) = — cot 62, 
give the interface shape outside the plates as: 
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while the contact angle conditions (0) — — cot 9i and 
h'2 (d) = cot 62 give the interface shape between the two 
plates to be: 
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Because of the interfacial deformation given by (2)- 
(4), the plates are now subjected to a capillary pressure 
that results in a horizontal force on the plates (there is 
no resultant horizontal surface tension force because its 
components on either side of a plate cancel exactly). This 
force may act either to bring them together or to pull 
them apart depending on the two contact angles 0i and 
62- The value of the horizontal force per unit length, Fh, 
can be calculated by integrating the hydrostatic pressure 
along each of the wetted sides of one of the plates (say 
the one on the left in fig. 2]) and taking the difference as 
follows: 
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so that we have: 
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2 V sinh2(d/L,) 
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where the sign convention is such that < Q corre- 
sponds to attraction between the plates. Typically this 
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FIG. 5: A typical force-separation curve for a wetting and 
non- wetting plate and a liquid-gas interface. Here, 9i — 2n/3 
and 02 = 7r/4 and we observe repulsion at large separations 
as well as attraction at short range. 



force is either attractive for all plate separations or re- 
pulsive at large separations and attractive at short sep- 
arations (with an unstable equilibrium at an interme- 
diate distance). An example of a force-displacement 
curve in the latter case is shown in fig. [5] Equation ([6]) 
can be used to show that repulsion is possible only if 
cot 9i cot 6*2 < 0, i.e. if one plate is wetting and the other 
non-wetting, as follows. First, we note that, provided 
that cot 01 7^ — cot 612, the fact that F/i — + as c? — > oo 
implies that repulsion can only occur if Fii has a maxi- 
mum value somewhere, since asd^O, ^X— >— oo. A 
simple calculation shows that the only turning point of 
the function /(^) = (cot 6*1 cosh^ + cot02)/ sinh^ occurs 
at ^ = ^* where cosh^* = — cot 02 / cot 0i , which only 
has a real solution £_* provided that cot 6*1 cot 02 < 0. 
When cot 0i = — cot 02 , the short range attraction does 
not exist and instead there is repulsion at all displace- 
ments, but the result that repulsion can only occur when 
cot 01 cot 02 < still stands. 

The result just derived shows that vertical plates at a 
liquid-gas interface will attract if they have like menisci 
and repel otherwise, as we saw in section HIl with floating 
objects. However, the physical mechanism here is subtly 
different. We can no longer argue in terms of one plate 
following the meniscus imposed by the other as these 
plates do not float, meaning that there is no analogue of 
the condition of vertical force balance in this scenario. 
Instead, we must consider the effects of hydrostatic pres- 
sure that result from the deformation of the interface, 
as detailed by Poynting and Thomson®. Here, we give 
an abbreviated version of their argument in terms of the 
configurations shown in fig. [5] (a) and (b) in which plates 
of like wettability are at the interface. In (a), the up- 
per portion of the central column of fluid is at a pressure 
lower than atmospheric pressure, Patm, because of the 
curvature of the interface and the two plates attract be- 
cause of the excess of atmospheric pressure. Similarly 
in (b), the pressure in the outer fluid is greater than 
Patm, again because of the sign of the interfacial cur- 
vature, and so there is an excess pressure (or capillary 
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FIG. 6: Typical interfacial profiles for the different configura- 
tions of plate wettabilities, (a) Two wetting plates, (b) Two 
non-wetting plates, (c) A wetting and non-wetting plate at 
intermediate separation, (d) A wetting and non-wetting plate 
at short range. 



suction) causing the two plates to attract. The situation 
is more complicated when one plate is wetting and the 
other non-wetting. At intermediate and large displace- 
ments, the interfacial displacement of the central column 
at the points where it touches a plate is smaller than it 
is on the other side of the same plate (as shown in fig. [S] 
(c)) due to the constraint that the central interface must 
pass through z — Q (rather than just being asymptotic 
to 0) to allow it to satisfy both contact angle conditions. 
From ([5]), this is enough to show that the two plates re- 
pel. When the two plates come close to contact, however, 
a relatively large change in the gradient of the intermedi- 
ate meniscus is required between the two plates, inducing 
a large curvature. Since the curvature of the interface is 
proportional to its height, this in turn means that the dis- 
placement of the interface must be large in this region, 
reversing the sign of the force and leading to attraction 
(see fig. m (d)). The exception to this is the case where 
cot 01 = — cot 02 , as noted in the discussion following ([6]) . 
Here, very little curvature is required since the contact 
angles are precisely complementary and so, in turn, very 
little displacement of the interface is necessary. Thus, 
the interfacial displacement of the inside of the plates 
is smaller than that on the outside of the plates and so, 
again using ([5]), we see mutual repulsion between the two 
plates, regardless of the displacement between them. 

The argument just given has often been invoked to ex- 
plain the Cheerios effect (see, for example, the answer 
to problem 3.100 of Walker's book^ or perform a Google 
search on 'Cheerios effect'). While the effects of hydro- 
static pressure imbalances caused by interfacial deforma- 
tion are certainly important, the calculation considered 
in this section cannot constitute a complete explanation 
since it neglects the crucial fact that Cheerios, and other 
fioating objects, must satisfy a vertical force balance to 
be able to float. As we shall see in the next section. 
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FIG. 7: Geometry of a sphere lying at a liquid-gas interface. 
The shaded area represents the weight of liquid equivalent to 
the buoyancy force due to hydrostatic pressure acting on the 
spherei^. 



this changes the physics fundamentaUy for smaU par- 
ticles, leading to the buoyancy mechanism outhned in 
section llll providing the dominant effect. Incorporating 
this additional equilibrium condition, however, requires 
a slightly more involved calculation, which we shall now 
consider. 



V. 



THE CASE OF FLOATING OBJECTS 



Having seen the manner in which the force between two 
interfacial objects can be calculated in a somewhat arti- 
ficial geometry, we are now in a position to consider the 
scenario that is of most interest to us here: interactions 
between objects that are floating at a liquid-gas inter- 
face. The major difference between this and what has 
gone before is that we must now take into account the 
vertical force equilibrium of the object. This is actually a 
significant enough complication that, even using the lin- 
earised approach of section llVl progress can only be made 
numerically. However, for sufficiently small particles we 
may assume in addition that the interfacial profiles gen- 
erated by one or more object floating at the interface are 
sufficiently small that they may be superposed and so 
make progress analytically. 

This assumption was tacitly made in our discussion 
of the attraction and repulsion of objects at interfaces 
in section [Hi but was first introduced by Nicolson^ who 
used it to calculate the interaction force between neigh- 
bouring bubbles in bubble rafts. It was then applied to 
the calculation of the force between floating particles at 
an interface by Chan et al.^ who considered some simple 
illustrative particle configurations such as two horizon- 
tal cylinders floating near one another. Here we give 
their argument applied to the problem of determining 
the interaction between two identical spheres floating at 
an interface, primarily motivated by the simplicity of an 
experimental realization thereof. 



In the spirit of the Nicolson approximation, we first 
neglect the presence of the second sphere and consider 
the vertical force balance on an isolated sphere to deter- 
mine an approximate expression for the interfacial slope 
at the contact point, z'^ (see fig.E]). This is equivalent to 
determining the value of (/)c, as defined in fig. [71 because 
of the geometrical relation (pc ~ tt — 9 + arctan z^, but 
the calculation of z'^ will prove useful to justify Nicolson's 
approximation and to compare the results with those of 
section IIVI 

For this sphere to remain at the interface, its weight, 
^irpsgB? must be balanced by the component of sur- 
face tension acting along the (circular) contact line 
and the buoyancy force because of the displaced bulk 
fluid. The first of these is easily seen to be given by 
27ri?sin0c7sin(arctanz^) — iTT^RsiiKpcz'^i^ + z'^)~^/'^^ 
while the second is just given by the weight of water 
that would occupy the area between the wetted region 
of the sphere and the undisturbed interface, which is 
shown as the hatched area in fig. [71 To understand phys- 
ically this generalisation of Archimedes principle, notice 
that the liquid has no knowledge of the geometry of 
the object that is at the interface outside of its wetted 
perimeter. The liquid must therefore produce an up- 
thrust equal to that it would provide to an object filling 
the entire hatched region, which we know from the usual 
Archimedes result is just the weight of displaced liquid 
that would fill this volume. (For elegant rigourous deriva- 
tions of this result, the reader might consult Keller— or 
Mansfield et al}^). This volume can be calculated by 
splitting it into a circular cylinder of radius R sin (j)c and 
height Zc and a spherical cap of height R{1 — cos(/)c) and 
base i? sin 0c- This gives the buoyancy force to be: 
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The balance of these vertical forces may now be written 
explicitly as: 
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By substituting (pc = tt — 6 + arctan and keeping 
only those terms linear in z'^ we obtain an expression 
for z'^ sin (jfc accurate to linear order in the Bond number, 
B ^ R^/Ll: 

(2D - I 1 1 , \ 

Zc sin 0c = S cos 61 + - cos^ 6* = BE (9) 

\ o D / 

where D = Psj P- (As a consistency check, observe that 
Zc = when Q = 7r/2 and D = 1/2, which we ex- 
pect, since in this case the Archimedes buoyancy alone is 



enough to balance the weight of the sphere without any 
interfacial deformation.) 

Equation ([9|) contains two dimensionless parameters, 
B and S, whose physical significance it is important to 
discuss. The Bond number, B = pgR^ /j, is the most im- 
portant dimensionless parameter in this system. Physi- 
cally, it gives a measure of the relative importance of 
the effects of gravity and surface tension: large B cor- 
responds to large particles or small surface tension co- 
efficient - in both cases surface tension is inconsequen- 
tial. The expression for the slope of the interface in the 
vicinity of the spherical particle given in ^ is valid for 
B ^ 1 (corresponding to a radius of ~ 1mm or smaller 
for a sphere at an air-water interface) in which case sur- 
face tension is very important. The other dimensionless 
parameter, S, can be thought of as a (non-dimensional) 
resultant weight of the particle once the Archimedes up- 
thrust has been subtracted out. This physical inter- 
pretation arises naturally from the vertical force bal- 
ance condition ([5]) and (O since the resultant weight 
of the object (in the linearised approximation) is simply 

To calculate the interaction energy using the Nicolson 
approximation, we must also calculate the interfacial dis- 
placement caused by an isolated floating sphere, which is 
determined by the hydrostatic balance 7V^/i — pgh - the 
co-ordinate invariant statement of equation ([T]). With 
the assumption of cylindrical symmetry, this becomes: 
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with the boundary conditions that h —^ as r —^ oo and 
h'{r = R sin ipc) = z'^. Equation (fTO|) has a solution in 
terms of modified Bessel functions of the first kind^^ and 
of order i, Ki{x): 



~z'^ sin (f)cRKo{r/Lc) 
(11) 

where we have used the asymptotic result (see, for ex- 
ample, Abramowitz and Stegun^^) that Ki{x) w l/x for 
X <C 1 to simplify the prefactor. 



B. Two interacting particles 

Having calculated the effective weight of a sphere at 
a deformed interface as 2tt^RB'S (with S as defined in 
equation ^) as well as the interfacial deformation caused 
by the presence of a single sphere we are now in a posi- 
tion to calculate the energy of interaction between two 
spheres. To leading order in B, this energy is the product 
of the resultant weight of one sphere and its vertical dis- 
placement caused by the presence of another sphere with 
its centre a horizontal distance / away. We may therefore 
write the energy, E{1), as: 
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and from this, the force of interaction is given by F{1) 



dl 
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Exercise: Carry out the above calculation for the 
case of two cylinders of infinite length lying horizontally 
and parallel to one another at an interface. To do this, 
first consider the interfacial profile caused by an isolated 
cylinder and show that this is given by 



z{x) = ~Lcz'^e-yiY>{~x/Lc) 



(14) 



when i? <C 1. Next use the linearised vertical force 
balance and the geometrical relationship (pc ~ tt — 9 + 
arctanz' to show that: 
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The resultant weight of the object can be found from 
the force balance to be « 2'^z'^. Use this and the inter- 
facial profile in to calculate the energy of interac- 
tion (per unit length) between two cylinders with centre- 
centre separation E(X). Show that: 
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with C defined as in ([T5 



It is important to emphasise that the calculation lead- 
ing to p2p relied on a number of assumptions. The first 
of these is that the particle be small enough that the 
total interfacial deformation is the sum of that due to in- 
dividual particles (by comparison with numerical results, 
Chan et al^ show that the expression derived by their 
method is essentially exact for Bond numbers B < 0.1). 
Furthermore, we have neglected the effect of capillary 
pressure acting on the particle to produce a horizontal 
force - calculated for the case of two vertical plates in 
section HVl It is not possible using the analysis presented 
here to include this effect since inherent in the Nicolson 
approximation is the assumption that the level of the 
interface is the same on either side of the particle. How- 
ever, the difference in interface heights on either side of 
the sphere only occurs at the next order in B so that 
using ((S]), we see that the contribution from capillary 
pressure also only enters at the next order in B without 
doing the detailed calculation. 

At large Bond numbers the analysis presented above 
breaks down raising the question: is it still the case that 
large particles interact because of their gravitational po- 
tential energy or does the capillary suction effect dis- 
cussed in section Hvl become more important? To answer 
this question conclusively requires the numerical solution 
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of the full problem since for large Bond number the inter- 
facial deflections are no longer small; a calculation that 
is beyond the scope of this work. Such a calculation has, 
however, been performed by AUain and Cloitre^'^, who 
showed that for _B ^ 1 the vertical displacement of two 
horizontal cylinders docs not change substantially as they 
move towards one another. This in turn means that the 
attractive force must result largely from pressure effects 
rather than the weight of the particles, with the crossover 
between these two regimes occurring for 10 < B < 100 
according to AUain and Cloitre's numerical results. The 
toroidal shape of a Cheerio complicates the notion of 
Bond number for a Cheerio but if we take an effective 
radius based on its volume R* = {R\R2Y^^ ~ 2.7mm 
(where Ri w 2mm and R2 ~ 5mm are the two radii of 
the torus) and Lc — 2.7mm for an air-water interface 
then _B w 1. This is within the regime where the grav- 
itational energy of the particles dominates the capillary 
suction due to the meniscus between them and so it is 
crucial that we account for the buoyancy effects to cor- 
rectly interpret the attractive force. 

Finally, we discuss briefly how the result in (fT3| fits to- 
gether with the physical interpretation of attraction and 
repulsion that we developed in section|TTl If they perform 
the exercise below, the reader will see that the sign of 
the interaction force between two non-identical spheres is 
governed by the sign of z'c^^ sin(/)i^' x z'}"^^ sin^c^'', where 
the superscript (i) labels the two particles. In particu- 
lar, the interaction is attractive if this product is positive 
and repulsive if it is negative. From pT|) . we see that 
this product has the same sign as the product of the gra- 
dients of the menisci in the neighbourhood of the two 
particles, and so we see that there is mutual attraction 
if the particles cause like menisci and repulsion if they 
have unlike menisci. This is precisely the picture that 
we saw in section |TT1 although we are now able to quan- 
tify the combination of contact angles and particle den- 
sities that give rise to the two different possibilities. We 
also note that the strength of the interaction decreases as 
the surface tension coefficient, 7, increases. This slightly 
counterintuitive result is actually a simple consequence 
of the fact that for higher values of 7 the deformation of 
the interface required to satisfy the vertical force balance 
condition for particle 1 is less and so the gravitational 
hill upon which particle 2 finds itself is smaller. 



Exercise: By repeating the analysis for two spheres 
but now allowing them to have different material prop- 
erties R, 9 and D, show that the sign of the interaction 
force for two (non-identical) spheres is determined by the 
sign of z'c^'^ sm(j)c^ x z'c'^'^ s\Q-(j3i \ where the superscript 
[i) labels the two particles (the interaction is attractive 
if this product is positive and repulsive if it is negative). 



VI. THE DYNAMICS OF FLOATING 
PARTICLES 

Having so far limited ourselves to calculating the force 
of interaction between particles, we shall now use this to 
answer some questions that arise naturally from observ- 
ing the motion of objects at an interface. Perhaps the 
most natural question to start with is "how fast do two 
spherical particles come together?" and it is this that we 
shall focus on in this section. 

To make the dynamical problem tractable, we assume 
that the vertical velocity of each particle as it moves along 
the meniscus is small enough that the vertical force bal- 
ance used to determine the horizontal force is satisfied; 
an assumption that is valid to leading order, as shown by 
Chan et al^ . 

We assume that the motion is overdamped so that a 
Stokes drag term* (for the viscous drag provided by the 
liquid) balances the attractive force between the two par- 
ticles given in p^ . This balance leads to the following 
equation of motion: 

&^T^J.Ra^ = -2^7i?i35/2i]2x^ ( M\ (17) 
at \Lc J 

where ^ is the dynamic viscosity of the liquid and a is 
a scaling factor to take into account the fact that the 
drag experienced by a particle at an interface is less than 
it would experience if completely immersed in the bulk 
fiuid. We expect that a « 1/2 although Danov et alw^ 
computed numerically the dependence of a on 9. Data 
obtained from observing two spherical particles with ra- 
dius 0.3 mm as they move under each other's influence 
at the interface between air and water is shown in flg. 
[H This data was collected from a time lapse movie (one 
frame per second) of the motion taken with a digital cam- 
corder, which was then analysed using image analysis 
software (Image J, NIH). Here the non-dimensional resul- 
tant weight of the particles, S, is difficult to measure 
experimentally because of its dependence on the contact 
angle 9. We thus appear to have two unknown param- 
eters in this model (a and S). Fortunately, since only 
the ratio Y? /a appears in ([TT]) we can fit the numerical 
solution of pT|) to the experimental data presented in fig. 
|S]via this one parameter. 

In this case, we are assuming that the particles are 
sufficiently small that their inertia may be neglected en- 
tirely. In fact, the approach we have adopted is better 
suited to such small particles since the Bond number in 
this situation is also very small and so the expression for 
the interaction force in is effectively exact. Also, 
if the typical distance between particles is small (com- 
pared to the capillary length, Lc) then we are again able 
to use the asymptotic formula Ki{x) for a; <i; 1 to 

approximate the modified Bessel function that occurs in 
the force-law (IT5|) . Substituting this into ([TT]) allows us 
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FIG. 8: Experimental data (points) compared to the solu- 
tions of the dynamical equation (17) (solid line) for two iden- 
tical spheres of radius 0.3mm interacting via flotation forces 
at an air-water interface. The theoretical curve is calculated 
by solving (17) numerically with 7 = 0.0728 Nm~\ = 0.001 
Pa s and Y? ja = 0.673 used as the fitting parameter. (S is 
the non-dimensional resultant weight of the spheres and a is 
the ratio of the drag that the particle feels at the interface to 
that which it would experience in an unbounded fluid). The 
dashed line gives the asymptotic result (|18p . 



to solve for Z(i), giving: 



3/iQ; 



(18) 



so that the time taken for two spheres to come into con- 
tact is given by: 



(19) 



In fig. 8, we see that the asymptotic form p8)) is a rea- 
sonable approximation over the last 2 seconds before con- 
tact. 

To verify our assumption that the motion is slow 
enough for us to be able to neglect the particle's accel- 
eration, we look at the ratio of particle inertia to Stokes 
drag. This ratio TZ is initially small and increases as the 
particles come closer, reaching a maximum when they 
come in contact. Using the asymptotic expression for 
Ki {x) , we have: 



Tic 



2 jpL.B'^/^j:^ 



27a2 



(20) 



For spherical particles of radius less than 0.3mm (or 
B = 0.01), TZ < 0.1 throughout the motion, which is 
sufficiently small that our approach is self-consistent. 



VII. AMPHIPHILIC STRIPS WITHOUT 
CHEMISTRY 

As a final illustration of the principles that we have 
applied to the problem of the attraction of interfacial ob- 
jects we shall consider briefly the equilibrium of a single 



FIG. 9: A single two-dimensional strip (of weight W per unit 
length in the direction perpendicular to the page) floating 
at an interface, with the hatched area indicating the area of 
displaced fluid, the weight of which is equal to the buoyancy 
force on the strip. The asymmetric equilibrium position here 
is possible because of an off-centre position for the strip's 
centre of mass. 



strip (of weight W per unit length) floating horizontally 
at an interface, as shown in fig. [HI This problem was 
studied extensively by Mansfield et alM- who also inves- 
tigated the interactions between two such strips. Here, 
however, we content ourselves with studying the simpler 
problem of a single strip but with a slight twist: we take 
the centre of mass line of the strip to be displaced from 
the strip's centre line by a distance (1 — (3)b. As we shall 
see, this is enough to break the symmetry of the prob- 
lem and thus allow the strip to float at an angle a to the 
horizontal. 

Because the strip is assumed to be infinitely thin, the 
concept of contact angle that we met earlier is not well- 
defined in this problem. Instead, the effective contact 
angles that the interface makes at the points of contact 
with the strip, </'i,2, and the angle a are determined from 
the condition that the strip be in equilibrium, which leads 
to three equations for the three unknowns. The first two 
of these are the conditions of horizontal and vertical force 
balance, which we have encountered previously and here 
may be written: 



= 7(cos(/)i - cos(/)2) + ^{zl ~ zl) 
for horizontal equilibrium and: 



W = 7 (sin 01 -f sin 02 ) — pgC^b cos a) 



Zl + Z2 



(21) 



(22) 



for vertical equilibrium - the second term on the right 
hand side in (j22p simply being the weight of liquid dis- 
placed by the strip, which is shown as the hatched area 
in fig. [HI In addition, there is the condition of torque 
balance, which we did not encounter for the equilibrium 
of the floating sphere (it is automatically satisfied for ob- 
jects of circular cross-section). Here, however, this con- 
dition is crucial as it gives a third equation by which 
to determine the three unknowns. By taking moments 
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about the centre of mass, we have: 



Finally, the torque balance condition ([23|) yields: 



= 76((2 - ^)sin((/)2 - a) - /3sin((/)i + a)) 

.(2-/3)fc 

— pg z{s)sds (23) 

J-l3b 

where s is an arc-length co-ordinate measured along the 
strip from the centre of mass. 

In principle, these equations can be solved even for 
interfacial deformations that are not small following the 
strategy outlined by Mansfield et al}^ . However, for sim- 
plicity, we proceed here in the limit where the three an- 
gles are small and, therefore, the interfacial deformations 
are also small. In this limit zi ~ —Lc(f>i,Z2 ~ —Lc4>2 so 
that the vertical relation (1^^ becomes: 



W 

7 



1 



b 



(24) 



Then, the linearized version of (pij) is automatically sat- 
isfied and so we make use of the geometrical relation 
Z2 — zi = 2&sina to give: 



(25) 



26/Lc 



(2 - /?)02 - - 2a 
^ (l-/3)(/302-h(2-/3)0i) 



L 

b Cj)2 



(26) 



Together, equations ((24|) - (l26l) constitute a system of 
three linear equations in three unknowns, which can be 
solved simply by inverting a 3 x 3 matrix . Of particular 
interest here is the possibility that, for some values of /3, 
we may have 4>i and 4>2 of opposite sign and so plotted 
in fig. [TO] is the dependence of (t)\4'2 on (3 showing that 
for sufficiently small /3, (fix and (f)2 are of opposite sign. 
We can understand this more fully by solving the system 
(IMl)-(Eni) to give: 



2 {3Ll + 6bU + 46^ - 3b{Lc + b)/3){3Ll - 29 + 3b{L^ + b)P) 
A{L, + bY{b-^ + ibL, + iLlY 

I 



(27) 



From this and the restriction that /3 < 1, it is a simple 
6i02 < when /3 < /3c ^fc^-si^ 



matter to show that 



3b(fc+Le) 

and hence that (^i and (j)2 are of opposite sign for suf- 
ficiently small P provided that b > (3/2) ^/^Lc. For the 
parameters in fig. [TOl /3c « 0.28, with (j)2 changing from 
negative to positive around (3 = (3c and 0i positive for all 
(3 G [0, 1]. Physically, this is what we should expect since 
for small /?, the offset of the centre of mass causes a large 
torque and hence rotation of the object, forcing the far 
end of the strip to be displaced above the equilibrium liq- 
uid level and hence the interface here must be elevated, 
corresponding to 4>2 < 0. As the offset decreases, so does 
the torque associated with it and the displacement of the 
far end of the strip is diminished until, for sufficiently 
small offsets, the far end lies below the equilibrium liq- 
uid level and 02 > 0. 

The significance of the qualitatively different interfa- 
cial shapes at either end of the strip that we see when 
(/)i02 < is that in such cases one edge is able to attract 
drawing pins while the other repels drawing pins and 
attracts bubbles! Particles or molecules with this type 
of behaviour are often termed amphiphiles and occur in 
detergents along with many other applications-^. How- 
ever, such amphiphilic particles are usually constructed 
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FIG. 10: The dependence of <^i(?i2 on the offset of the strip's 
centre of mass, given by (1— /3)fe, for /3 G [0, 1] (the graph being 
symmetric about f3 — 1). Here the strip has non-dimensional 
weight per unit length W/'y = 0.1 and half- width b/L^ — 2. 



by treating the two edges chemically to induce the very 
different behaviour: here, we have shown that this may 
also be achieved by altering the density of the strip so 
that the centre of mass of the object is displaced. 

As well as the possible industrial applications that par- 
ticles such as this could have, this physical amphiphile 
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FIG. 11: Example of a simple amphiphilic strip, (a) Simple 
design of the amphiphilic strip mentioned in the text, (b) 
Photograph of such an amphihilic strip bonding a drawing 
pin (left) and the cap of another drawing pin (right), which 
would normally be mutually repulsive. Using such an am- 
phiphile allows particles that would otherwise experience mu- 
tual repulsion to be maintained stably at an almost arbitrary 
separation. 



lends itself to a much simpler classroom demonstration 
than might be possible with chemical amphiphiles. Using 
a piece of sticky tape doubled back on itself to form the 
strip, see fig. [11] (a), and a short length of wire (part of 
a paperclip for example) inserted between the two sides 
of tape, it is possible to make a strip with edges that 
deform the interface in manifestly different ways. An 
experimental realization of this is shown in fig. [TT] (b), 
which demonstrates that such particles are an effective 
means by which particles that would otherwise be mu- 
tually repulsive can be coaxed into maintaining a finite 
equilibrium separation. 



VIII. DISCUSSION 

In this article we have investigated an aspect of every- 
day life that may previously have escaped our readers' 
notice - the propensity of fioating objects to aggregate. 
Having couched the mechanism for this effect in terms 
of the simple physics of particles trapped at a deformed 
interface feeling the effects of their weight (or buoyancy), 
we have shown how approximate methods can lead to 
quantitative descriptions of the magnitude and dynam- 
ical nature of the interaction. We now conclude with a 
brief discussion of some of the scenarios in which this ef- 
fect has found application and point out some potential 
research directions. 

As is often the case, the first place to look beyond 
the kitchen and the lab is the natural world. There are 
likely many instances where a variant of the Cheerios ef- 
fect is used by some species or another. Here, however, 
we choose to highlight one particular application: wa- 
ter walking creatures that can be found on the surface 
of many ponds. Many of these animals rely on surface 




FIG. 12: The interfacial deformation caused by a water-spider 
{Dolomedes triton). Here, the interface is visibly depressed 
by the spider's weight acting on each leg. Image courtesy of 
Robert B. Suter, Vassar College. 



tension to prevent them from drowning as their weight 
can be supported by interfacial deformations (as is shown 
in fig. fTS]) . However, when they try to climb out of the 
pond, they become reluctant victims of the Cheerios ef- 
fect since this generally requires climbing up the meniscus 
- against gravity. Recent observations by Hu and Bush- 
suggest that some insects, such as Mniovelia Kuscheli, 
may overcome this difficulty by pulling up on the inter- 
face with their front legs and pushing down on it with 
their hind legs, effectively shifting their center of gravity. 
They thus become the natural exemplar of the mechan- 
ical amphiphile discussed in the last section as they are 
now attracted to the wall via the Cheerios effect. 

There remain many interesting questions that we have 
not been able to answer in this article, many of which 
are certainly amenable to investigation in the classroom 
or the laboratory. For example, it would also be worth- 
while to better understand the way in which two differ- 
ent types of particles that are mutually repulsive interact. 
With a mixture of light and heavy particles, for example, 
clusters of like density particles form - all things (other 
than the particle density) being equal. This segregation 
phenomenon would undoubtedly find many applications 
within science and industry. Such a study could then 
naturally be extended by the inclusion of amphiphilic 
particles, which would allow the user to dictate a finite 
equilibrium spacing between two otherwise repulsive ob- 
jects. Another example is afforded by hair on water: a 
flexible hair floating parallel to a planar wall might be 
expected to bend as it is attracted to the wall since the 
attractive force would be greatest for those parts of the 
hair closest to the wall and much less signiflcant for those 
further away. This motion has yet to be studied from ei- 
ther an experimental or theoretical point of view. The 
inspiration from the kitchen, industry and nature seems 
almost overwhelming. 
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